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LONG-TERM  GOALS 

Smoothed  Particle  Hydrodynamics  (SPH)  is  a  novel  meshless  numerical  method  that  is  being 
developed  for  the  study  of  nearshore  waves  and  Navy  needs.  The  Lagrangian  nature  of  SPH  allows 
the  modeling  of  wave  breaking,  splash-up,  and  the  subsequent  fluid  turbulence,  which  in  large  part  is 
comprised  on  coherent  turbulent  structures.  Including  sediment  transport  allows  the  methodology  to 
be  applied  realistically  within  the  surf  zone. 

OBJECTIVES 

The  objectives  of  this  project  are  to  improve  the  SPH  model  for  use  in  unraveling  the  physics  of 
breaking  waves,  including  the  description  of  the  wave-induced  turbulence  and  sediment  transport 
within  the  surf  zone.  Coupling  the  SPH  model  to  a  vertically-averaged  more  computationally  efficient 
model  for  larger  areas  would  provide  a  computationally  useful  model  of  the  full  surf  zone. 

APPROACH 

The  approach  is  based  on  improving  various  aspects  of  the  JHUSPH  code;  applying  the  code  to  more 
validation  tests;  and  to  examine  in  some  detail  new  aspects  of  the  model  by  applying  it  to  different 
situations.  Explorations  of  hydrid  models,  that  is,  coupling  the  SPH  particle  model  to  conventional 
finite  difference  models,  such  as  the  Boussinesq  model,  FUNWAVE,  are  being  explored. 

WORK  COMPLETED 

FY06 

•  Modifications  to  the  JHUSPH  code  this  calendar  year  involve  implementing  a  new  time 
integration  scheme,  based  on  Beeman’s  (1976)  method,  which  is  0(At4),  which  is  two  order  higher 
than  the  previous  method,  based  on  the  Verlet  (1967)  algorithm. 

•  Validation  of  the  SPH  hydrodynamics  calculation  has  been  carried  through  comparisons  to 
laboratory  experiments. 

•  Development  of  a  Large  Eddy  Simulation  equivalent  for  the  turbulent  contributions  to  the  sub¬ 
particle  suspended  sediment  transport. 
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•  Validation  of  the  model  for  suspended  sediment  transport  for  laboratory  experiments  of 
oscillatory  flow  over  ripples. 

•  Implementation  of  a  scheme  to  identify  coherent  turbulent  structures  in  breaking  waves. 


RESULTS 


Code  Improvements:  The  code  was  improved  this  year  with  the  implementation  of  a  variety  of  new 
algorithms,  such  as  a  better  time-stepping  method  (based  on  Beeman,  1976)  for  a  higher  accuracy;  sub¬ 
particle  scaling  for  turbulent  suspended  sediment  transport  to  account  for  turbulence  at  length  scales 
smaller  than  the  size  of  the  SPH  particle  (consistent  with  the  sub-particle  scaling  used  for  the 
hydrodynamics  turbulent  closure  scheme);  and  included  a  new  kernel  (weighting  function),  taken  from 
radial  basis  functions  (Wendland,  2005),  which  requires  much  less  computation  than  the  cubic  spline 
function  we  have  used  earlier.  In  addition,  the  decision  has  been  made  to  implement  a  more  robust 
numerical  approximation  into  the  SPH  method.  Traditional  SPH  uses  a  weighting  function  based  on 
dispersed  fluid  “particles”  or  nodal  points  to  interpolate  data.  It  has  been  shown  that  this  approach  is 
similar  to  approximating  a  tilted  plane  with  a  constant  value.  Moving  Least  Squares  (MLS)  has  been 
suggested  by  Dilts  (1999)  and  others  as  a  more  exact  approach  to  the  approximation,  which  allows  a 
local  polynomial  approximation  to  the  plane.  This  is  being  implemented. 

Hydrodynamics  Validation:  SPH  is  a  robust  model  that  provides  realistic  results  for  nearly  all 
examples  of  free  surface  flows.  Despite  that,  very  few  details  comparisons  to  hydrodynamics  are 
extant  in  the  literature;  in  part,  this  is  due  to  the  need  to  have  high  resolution  for  these  comparisons. 

We  have  been  carrying  out  validation  tests  here  for  waves  on  beaches  (Synolakis,  1986  and  Ting  and 
Kirby,  1995)  and,  with  colleagues  at  the  University  of  Vigo,  Spain,  for  cases  of  a  dam  break  problem 
with  standing  water  in  front  of  the  dam. 

In  Figure  1,  a  comparison  is  made  between  JHUSPH  results  and  the  data  (given  by  blue  stars)  obtained 
by  Synolakis  (1986)  in  a  wave  tank  study  of  solitary  waves  on  a  beach.  The  motion  of  the  wave 
maker,  the  propagation  of  the  wave  across  a  horizontal  section,  and  run-up  of  the  wave  on  slope,  are  all 
accurately  modeled. 
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Figure  1  SPH  model  (in  green)  compared  to  data  from  Synolakis  (1986)  for  a 

shoaling  solitary  wave. 


Wave-breaking  turbulence:  A  plunging  breaker  on  a  beach  creates  a  tremendous  amount  of 
vorticity.  It  was  observed  in  our  previous  numerical  studies  that  the  vertical  vorticity  appears  in  the 
breaking  region  concurrently  with  the  breaking  event  and  persists  after  breaking,  evolving  into 
coherent  structures,  such  as  obliquely  descending  eddies. 

To  identify  these  coherent  turbulent  structures,  a  method  used  in  turbulence  is  adopted.  The  velocity 
in  the  small  neighbourhood  of  each  particle  is  assumed  to  vary  linearly,  namely, 

L'i  —  AijXj 

where  Aij  are  local  components  of  strain  rate  and  xj  are  local  coordinates.  The  eigenvalues  of  the 
tensor  A  can  be  obtained  from  the  characteristic  equation: 

A:s  +  PX2  ■+  QX  +  R  -  0 
P  -  -At, i  Q  ~  and  tt~=  A\ 

P,  Q,  and  R  are  the 
invariants  of  the  tensor  A. 

While  vorticity  can  be 
precisely  defined  as  the  curl  of 
vorticity,  there  is  no  exact 
definition  of  a  vortical 
structure.  It  has  been  shown 
(cf .  Bernard  and  Wallace, 

2002)  that  in  regions  of 
positive  values  of 
discriminant  D  of  the 
characteristic  equation  above, 
the  local  streamlines  indicate  a 
swirling  motion  about  the 
point.  Mathematically  the 
matrix  A  will  have  one  real 
and  two  complex  eigenvalues 
when  D  is  positive.  Hence 
the  magnitude  of  the 
imaginary  part  of  the 
eigenvalue  is  used  as  a 
criterion  to  look  for  vortical 
structures. 

For  each  water  particle  in  the  simulation  domain,  the  eigenvalues  of  A  are  computed.  Iso-surfaces  of 
the  magnitude  of  the  complex  eigenvalue  equal  to  some  small  positive  number  are  then  plotted  to  look 
for  vortices.  Figure  2  shows  an  iso-surface  visualization  of  vortices  corresponding  to  the  touchdown 
and  rebound  of  a  plunging  breaker  jet.  While  all  the  vortical  structures  shown  in  the  figure  above  are 
plotted  in  gray,  some  are  rotating  in  one  direction  and  their  neighbors  in  the  opposite  direction.  These 
vortical  structures  show  considerable  alongshore  (Y)  consistency,  yet  there  is  clearly  3-D  structure 
evolving,  turning  into  obliquely  descending  eddies. 


Figure  2  Coherent  vortical  structures  under  breaking  wave.  X 
increases  towards  the  beach,  Y  is  along  shore  distance,  and  red 
denotes  the  bottom. 


Suspended  Sediment  Calibrations 


Our  earlier  examinations  of  suspended  sediment  over  ripples  (Zou  and  Dalrymple,  2004,  2006)  has 
been  extended  by  the  development  of  a  parallel  code  version  of  oscillatory  flow  in  a  U-tube  (actually 
modeled  by  fluid  in  a  horizontal  tube  bounded  by  moving  pistons).  A  sub-particle  scheme  (like  LES) 
was  derived  for  the  turbulent  diffusion  term  to  resolve  the  suspended  sediment  flux  in  eddies  that  are 
smaller  than  the  particle  “size”  (the  size  of  the  compact  support  for  the  SPH  weighting  function). 

In  Figure  3,  the  comparison  of  the  JHUSPH  computed  vertical  distribution  of  volumetric  sediment 
concentration  over  the  crest  and  trough  of  wave-induced  ripples  of  length  X.  The  model  results 
compare  satisfactorily  with  the  experiments  (Ribberink  and  Al-Salem,  1989). 
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Figure  3  Volumetric  concentation  over  ripple  trough  and  crest  as  a  function  of  dimensionless 
elevation  above  the  bed.  Dimensionless  elevation  is  elevation  divided  by  ripple  length. 


International  Collaborations 

The  JHUSPH  project  at  Johns  Hopkins  University,  which  we  consider  as  an  open  source  effort,  has 
benefitted  from  collaborations  with  indiviuals  from  a  number  of  European  universities.  This  year,  we 
have  hosted  Dr.  Andrea  Panizzo  from  La  Sapienza,  Rome;  Giacomo  Viccione  from  University  of 
Salerno,  Italy;  and  Dr.  Moncho  Gomez-Gesteira  from  the  University  of  Vigo,  Spain  for  research  on 
SPH.  Two  of  these  individuals,  plus  former  Post-Doc  and  other  visitors,  have  played  a  major  role  in 
the  development  of  a  Special  Interest  Group  of  the  European  Research  Community  on  Flow, 
Turbulence,  and  Combustion  (ERCOFTAC)  called  SPHeric  (http://cfd.me.umist.ac.uk/sph/).  In  fact, 
over  20%  of  the  steering  committee  of  this  group  have  worked  on  SPH  at  Hopkins,  and  Johns  Hopkins 
University  (through  the  P.I. )  is  a  foreign  member  of  SPHeric. 


IMPACT/APPLICATIONS 


Smoothed  Particle  Hydrodynamics  is  proving  to  be  a  competent  modeling  scheme  for  free  surface 
flows  in  two  and  three  dimensions.  Coupled  with  another  wider-area  wave  model,  such  as 
Boussinesq,  a  hybrid  SPH  model  would  provide  a  large,  highly  resolved,  look  at  an  entire  surf  zone. 
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